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We discuss the advantages of using the momentum space lattice method for studying the phase diagram of 
quenched QED4. Preliminary results of a numerical simulation are presented. They indicate that the method 
avoids the contamination by 4-fermi interactions which plagues the conventional position space non-compact 
formulation. 



At last year's conference, we introduced the 
momentum space lattice method and discussed its 
general features |Q . Here we would like to report 
on its application to quenched QED4. Before go- 
ing into the details of the simulation, however, it 
is worth explaining the advantages of the method 
for the problem at hand. We will therefore start 
with a brief reminder of earlier analytical and nu- 
merical results and motivate from there the intro- 
duction of the momentum space idea. 

The analytical work of Miransky, Bardeen, Le- 
ung and Love (MBLL) indicates that the phase 
diagram of quenched QED should be described in 
the {a,G) plane, where a is the usual QED cou- 
pling constant and G is the strength of a 4-fermi 
interaction. From the ladder Schwinger-Dyson 
equation, it is in fact possible to determine the 
shape of the critical line and to estimate the value 
of the critical exponents along this line They 
vary continuously from their mean-field value in 
the pure 4-fermi theory to an essential singular- 
ity at the MBLL point. Numerical simulations 
of quenched non-compact QED4 appear to pick- 
up a particular point along this line {a = 0.44q;c, 
G = 3.06), implying that non-compact QED ac- 
tually has to be interpreted as a "mixture" of 
QED -|- 4-fermi. The values of the critical ex- 
ponents at this point {5 = 2.2(1), /3 = 0.78(8), 
v = 0.68(2)) were confirmed in more extensive 
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simulations Q and recent results strengthen the 
hypothesis of power-law scaling ||]. All of this 
however raises an important question: can we 
find a simulation technique that would allow us to 
move along the critical line and expose the physics 
of "pure" QED? We certainly expect that if such 
a technique can be found in the quenched case, it 
will play an important role in future simulations 
of the full theory. 

A key point for carrying-out the above program 
would be to identify precisely the origin of 4-fermi 
interactions in non-compact QED with action: 



(1) 



An interesting candidate is revealed by using the 
quadratic nature of Sq and computing exactly 
the integration over the gauge fields. This opera- 
tion leads to an effective action of the form: 

■^eff. = j'Y^{i'i')x{i'i')x+,, 



exp[-l/(0)/2/3]5i.(C/ = 1) 



where we have assumed 
D^,v{x -y) = 5f,i,V{x - y) 



(2) 



(3) 



The first term represents a 4-fermi interaction, 
the second is a tadpole improved kinetic energy 
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term, and the dots stand for the current-current 
interaction as well as 3-body and higher order in- 
teractions. The form of (2) suggests that in order 
to get an effective action which is closer to "pure" 
QED, we should try to get rid of the 4-fermi term. 
In fact, recent results on quenched models seem 
to indicate that this term, when present, is very 
efficient in masking the rest of the dynamics: crit- 
ical exponents appear to be largely insensitive to 
the choice of the potential V There is no 
objection in principle to subtracting the 4-fermi 
term directly from the original action. It is gauge 
invariant and has all the appropriate symmetries 
(in fact, it is just the strong couphng limit of (1)). 
However, in practice, this turns out to be rather 
difficult. One could think for example of subtract- 
ing the 4-fermi interaction (multiplied by A^) by 
adding to the action: 

where 0^ is a new vector field which is freely inte- 
grated over between and 27r. Notice the relative 
-I- sign in (4) opposite to (1). This technique is 
however very slow numerically. In many respects, 
the problems are analogous to those encountered 
in simulations at finite chemical potential, where 
a term similar to (4) would appear at first order 
in an expansion of the lagrangian in powers of fi 
(in this analogy, A is identified with the chemical 
potential /x). As in the chemical potential case, 
only small values of A can be considered here, so 
that the 4-fermi interaction could only be sub- 
tracted partially. A more radical modification of 
the model therefore seems necessary. 

An important aspect of the 4-fermi interaction 
in S^g is that it is directly related to the com- 
pact nature of the coupling of the fermions to the 
gauge field in the original model (as can be seen 
through the steps leading to (2)). On the other 
hand, the momentum space lattice method 0] has 
a non- compact coupling and therefore avoids the 
generation of 4-fermi interactions. In fact, the 
method simply uses the continuum lagrangian of 
QED. The model is investigated on a torus of 
length L and the momentum expansion is trun- 
cated. After gauge fixing (to the Feynman gauge) 



and Fourier expansion, the action reads: 

+ 5]x(fc)[7^(27rfc^) + p]x(fc) (5) 

k 
k,k' 

where, instead of the usual fields tJj and A^, we 
have used the following dimensionless variables: 

X{k) = L^/V(fc) 0^.{q) = eL^A^iq) (6) 

and /? — 1/e^ , p = mL. In the above formula, 
the fc^'s are integers or half-integers depending on 
whether the boundary conditions for the fermions 
are periodic or antiperiodic. The q^'s are integers 
and take all the possible values of the momentum 
transfer between k and k' . The gauge fields there- 
fore live on a lattice which is twice as big (in each 
direction) as the fermionic lattice. The method 
has several advantages among which is the fact 
that it docs not suffer from the fermion doubling 
problem. From the numerical point of view, the 
procedure is however more involved than conven- 
tional position space lattice methods, the rea- 
son being the non-locality of the coupling of the 
fermions to the gauge field. This is partially taken 
care of by repeated use of a fast Fourier transform 
which reduces the computer time requirements to 
0(A^log A^). When quoting our results below, we 
give the size of the "photon lattice" since it is the 
one used in the FFT and therefore determines the 
computer time requirements. This size in turn is 
constrained by the specifics of the FFT routine 
and the choice of the boundary conditions for the 
fermions. Up to now, we have considered lattices 
of size (B)**, (10)^ and (15)^. We have measured 
the mass gap and < ipip > by respectively us- 
ing point and volume sources on the momentum 
space lattice (both methods make use of the fact 
that the propagator should be diagonal in mo- 
mentum after statistical average). Identifying the 
mass gap M with the inverse of the correlation 
length, we expect at the critical coupling: 

M oc m^/f^^ (7) 

From our (10)^ measurements , we obtain the es- 
timates /3c = 0.8 and (3S/i^ ~ 2. We can go one 
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step further and write the equation of state in the 
form: 

mSa: = ^(]^) A/? = /?-/3. (8) 

Fixing (35 /v = 2. and varying (3c and v, we obtain 
a best fit to our fuU data sample for (3c — 0.79 
and V — 1.22. This value of v is anomalously 
large and indicates that we are indeed probing 
the physics of a non-trivial theory (By compari- 
son, the value of v obtained on a position space 
lattice was v — 0.68 and it would he v — 0.5 in 
mean field theory). Similarly, the value (35 /v "^2 
is consistent with Miransky scaling. It is impor- 
tant to complete this picture by independent mea- 
surements of the critical indices (3 and 5. We mea- 
sured < ■i/'V' > on lattices of size (10)'' and (15)^ 
and found respectively (5 ~ 2.0 and 5 ~ 1.7. The 
trend here is interesting and indicates that as the 
size of the lattice is increased, it becomes possible 
to detect exponents which are further away from 
their mean- field value (see fig. 1). Since, on a 
position space lattice, distinguishing between the 
critical behaviour of QED and a Nambu Jona- 
Lasinio model is a difficult task, it is interesting 
to see how the momentum space lattice behaves 
in this respect. For that purpose, we ran simu- 
lations in which we replaced [InqY in equation 
(5) by a constant, so that the long range QED 
potential becomes a point-like interaction. The 
estimates for 5 are now 5 = 3.5 ± 0.2 on a (lO)'' 
lattice and 5 = 3.2 ±0.1 on a (15)4 lattice. These 
can very easily be discriminated from the QED 
values given above and we reach the conclusion 
that on a momentum space lattice, the shape of 
the potential really matters (at the difference of 
the position space situation [^). 

The results presented above are obviously 
very preliminary and more detailed analyses are 
clearly needed. However, they all indicate that 
the momentum space lattice method is successful 
in removing 4-fermi interactions from QED sim- 
ulations. The trends that we have established as 
a function of the lattice size, open the way for 
more detailed simulations of QED on larger lat- 
tices and the extraction of "infinite volume" crit- 
ical exponents. The method of course also has its 
limitations: it is computationally very expensive 
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Figure 1. Chiral condensate vs bare quark mass 
at various values of the coupling in the vicinity 
of the transition on a (10)^ lattice (lower group) 
and on a (15)** lattice (upper group) 



due to the non-locality of the interaction and a 
careful analysis of its gauge invariance properties 
is needed. It remains however, so far, the only 
method we have with the potential to map out 
the entire critical line of the quenched model. 
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